IL NUOVO CIMENTO 



Vol. ?, N. ? 



? 



GRB progenitors at low metallicities 

R. HlRSCHIlfll), G. MEYNETlpI), A. MAEDER(|21) 

Physics and Astronomy, University of Basel, Klingelberstr. 82, 4056 Basel, Switzerland 
Geneva Observatory, Ch. des Maillettes 51, 1290 Sauverny, Switzerland 



Summary. — We calculated pre-supernova evolution models of single rotating 
massive stars. These models reproduce observations during the early stages of the 
evolution very well, in particular Wolf-Rayet (WR) populations and ratio between 
type II and type Ib,c supernovae at different metallicities (Z). Using these models 
we found the following results concerning long and soft gamma-ray burst (GRB) 
progenitors: 

-GRBs coming from WO-type (SNIc) WR stars are only produced at low Z 
(Zlmc or lower). 

-The upper metallicity limit for GRBs is reduced to Z ^ 0.004 (SMC) when 
the effects of magnetic fields are included. 

-GRBs are predicted from the second (and probably the first) stellar generation 
onwards. 

PACS 00.00 - By the way, which PACS is it, the 00.00? GOK.. 
PACS — . .... 



1. Introduction 

Gamma-Ray bursts (GRBs) can be divided into two main types: a) short and liard 
and b) long and soft (see lT^ . Tiie long soft GRBs have now been firmly connected to 
supernova (SN) explosions of the type Ic (GRB 980425-SN 98bw, GRB 030329-SN 03dh, 
GRB 031203-SN 031w, GRB 060218-SN 06aj; see for example IH El). The leading model 
for GRB-SN explosions is the coUapsar model this model, the progenitor star 

must first lose its hydrogen rich envelope during its lifetime but retain enough angular 
momentum in its core. The star must also form a black hole when it dies. The large 
amount of angular momentum induces the formation of an accretion disk around the black 
hole. Then accretion from the disk onto the black hole powers bipolar jets. These jets 
can only pierce through the outer part of the star and be detected if the star is compact. 
This explain why the progenitor star has to lose its hydrogen rich envelope during its 
lifetime. Additional constraints come from the fact that close-by SN associated GRBs 
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are found (|a;|l6|lll) in low mctallicity {Z) environments, corresponding to the metallicity 
of the Magellanic Clouds {Z ^ 0.004 — 0.008). Low metallicity environments are also 
favoured from the theoretical point of view because lower metallicity means weaker mass 
loss and therefore smaller angular momentum loss (0). 

Several studies were devoted to finding which stars can lead to GRBs |3 S H IS ll^ 
I3:ll^l20[l. The key ingredients of the progenitor models in the context of GRBs are the 
treatment of rotation, magnetic fields, binarity and metallicity. In this work, we look 
at GRB progenitors coming from single stars with models, which include the treatment 
of differential rotation and its transport in the interior. We focus on the metallicity 
dependence of the GRB rates and discuss the impact of magnetic fields on our results. 

2. Computer model & calculations 

The stellar evolution code used for this work is described in (0). Convective stability 
is determined by the Schwarzschild criterion. Overshooting is only considered for H- 
and He-burning cores with an overshooting parameter, aovcr, of 0.1 Hp. The mass loss 
is proportional to the square root of the metallicity. It is also dependent on the surface 
rotation velocity. The centrifugal force is included in the structure equations. Last but 
not least, the processes taken into account, which induce transport and mixing of angular 
momentum and matter, are meridional circulation and dynamical and secular shears. 
Magnetic fields are included only in the models discussed in the section on magnetic 
fields. For these models, the treatment of magnetic fields ifTsj) is included according to 
l(l0|) . Models were calculated at various metaUicities: Z = 0.04 (GC), 0.02 (Zq), 0.008 
(LMC), 0.004 (SMC), and 10"* (second stellar generation) (seeUll for more details). 

3. Theoretical GRB rates around solar metaUicities 

Our models predict the production of fast rotating BHs. Models with large initial 
masses or high metallicity end with less angular momentum in their central remnant 
with respect to the break-up limit for the remnant. Many WR star models satisfy the 
three main criteria (BH formation, loss of hydrogen-rich envelope and enough angular 
momentum to form an accretion disk around the black hole) for long and soft gamma- 
ray bursts (GRB) production via the coUapsar model (T^. Considering all types of 
WR stars as GRB progenitors, we predict too many GRBs compared to observation. If 
we consider only WO-type WR stars (type Ic supernovae as is the case for the GRB 
related SNe discovered so far, see the list in the introduction) as GRB progenitors, the 
GRB production rates are in much better agreement with observations (see Table Pi . 
In particular, WO stars are found only at low metaUicities (Z<Zlmc) ^md have masses, 
M > 50 Mq. 

4. Do the first stars produce GRBs? 

We calculated a series of models ai Z = 10^^ (which corresponds to the second stellar 
generation). In this series, the 85 Mq model becomes a WO-type WR star (see Fig. Q] 
left). The core of the 85 Mq model retains enough angular momentum to produce a 
GRB via the coUapsar model (see Fig. fright). SNe of type Ib,c and GRBs are therefore 
predicted to ensue from the death of single massive stars at very low metaUicities. 
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Table I. - Predicted rates [yr ^] of long soft GRBs coming from WO-type WR stars, as well 
as limiting masses at various metallicities . 
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5.99E-04 










Observed rate: 7^B|| = 3 10"'' - 6 10"" (depending 
Reference total SN rate: 7^sN ^ =: 


on the beaming 
; 710"^ (14) 


angle of GRBs) 



5. Models with rotation and magnetic fields 

Petrovic et al |^j) and Fryer and Heger Q) studied the evolution of single and binary 
massive stars with and without magnetic fields (dynamo theory from TBI). They find that 
both single and binary models without magnetic field can produce GRBs and both single 
and binary models with magnetic field considered in their studies cannot produce GRBs. 
Models with magnetic fields nevertheless better reproduce the observed pulsar periods. 

The most recent studies found a way to produce GRBs in models including magnetic 
fields. Indeed, 12(1) and iflfli) find that if a star initially rotates fast enough, it undergoes 
quasi chemical homogeneous evolution. The key idea is that WR stars are produced by 




Fig. 1. - Left: profiles of the chemical abundances at the pre-SN stage for the 85 Mq model 
sX, Z — 10"*. Right: profile of the specific angular momentum (solid line) for the same model. 
The red dotted line shows the minimum angular momentum necessary in order to form an 
accretion disk around a rotating black hole. The blue short dashed and green long dashed lines 
show the minimum angular momentum necessary for matter around a maximally rotating and 
a non-rotating black hole, respectively. 
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mixing and not by mass loss. Since WR mass loss rates are too high at metallicities 
around solar, GRBs can only be produced via this channel in low metallicity environ- 
ments, find that GRBs can be produced for Z < 0.004. This is somewhat lower 
than the observed metallicity of close- by GRBs. 

We calculated 40 Mq models a.t Z — Z^uc and Z = 10^^ with magnetic fields. At 
Z = ^SMC) we obtain similar results as lj20f) . At Z = 10~^, we found that 40 Mq 
models undergo quasi chemical homogeneous evolution (if the initial velocity is larger 
than 600 km s"^). With or without magnetic fields, GRBs are therefore expected to follow 
the death of massive stars from the second stellar generation onwards (and probably from 
the first). 

The treatment of anisotropic mass loss and the recent downward revision of solar 
metallicity are two important physical ingredients that have yet to be included in the 
models and that will help GRB progenitors retain more angular momentum. 
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Appendix A. 
Questions cind answers: 

Question by A. J. van Marie: In several GRB afterglows, silicon is seen in absorption 
at a velocity that indicates its presence in the progenitor wind. Is this possible for low 
metallicity stars? 

Ansiuer. Yes, pristine silicon is present in the stellar wind. Nearby GRBs probably do 
not have very low metallicity progenitors, in which case pristine silicon could be detected. 

Question by D. Fox: You have investigated models of isolated massive stars, but many 
massive stars are born in close binaries. Does this affect your results? 
Answer. The presence of a binary companion may cither increase mass loss or induce 
accretion (at different times) and spin-up or spin-down a star. This makes it difficult to 
estimate the final effect for GRBs and obtain quantitative predictions. 

Question by C. Thompson: Could you summarize how your model of angular momentum 
transport is constrained by observations of isotope enhancements in evolved stars? 
Answer. One important test is comparing the surface enrichment, for example in nitro- 
gen, between models and observations as a function of the position in the Hertzsprung- 
Russell diagram. 

Comment by T. Piran: According to your low metallicity paradigm, we need new es- 
timates of GRB rates based on the SFR for low metallicity stars. 

Answer. I agree with this comment. Langer, Norman and collaborators have started to 
study this point. 



